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Abstract 

Flutter of an elastic body of revolution spinning about its axis of symmetry is prohibited in the subcritical 
spinning speed range by the Krein theorem for the Hamiltonian perturbations. Indefinite damping creates 
conical domains of the subcritical fiutter (subcritical parametric resonance) bifurcating into the pockets of 
two Whitney's umbrellas when non-conservative positional forces are additionally taken into account. This 
explains why in contrast to the common intuition, but in agreement with experience, symmetry-breaking 
stiffness variation can promote subcritical friction-induced oscillations of the rotor rather than inhibit them. 
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1 Introduction 

Consider an autonomous linear gyroscopic system describing small oscillations in the discretized models of 
rotating elastic bodies of revolution considered in a stationary frame [TJ [51 [31 [H [51 [5] 

x + 2f7Gi+(P + Sl2G2)x = 0, x = M2"^ (1) 

where is the speed of rotation, P=diag(ci;^, w^, c^l, . . . , w^, w^J = P-'" is the matrix of potential forces, and 
G=blockdiag(J, 23, . . . , nJ) = — G-^ is the matrix of gyroscopic forces with 
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Due to rotational symmetry of the rotor and periodic boundary conditions the eigenvalues < < < • • • < 
< of the matrix P arc double semi-simple, that is each eigenvalue has two linearly independent 
eigenvectors [71 |H] . The distribution of the doublets as a function of s is usually different for various bodies 
of revolution. For example, LOg = s corresponds to the spectrum of a circular string [6j. Nevertheless, there can 
exist isospectral bodies, because "one cannot hear the shape of a drum" [9]. 

Separating time by the substitution x = uexp(Ai), we arrive at the eigenvalue problem for the operator Lq 

Lo(r2)u (lA^ + 2f7GA + P + n^G^)u ^ 0. (3) 

As a consequence of the block-diagonal structure of the sparse matrices G and P, the eigenvalues of Lq are 

A^ = iuig ± isft, AJ = —iuJs T *sr2, (4) 

where bar over a symbol denotes complex conjugate. Rotation causes the doublet modes ±iLUs to split. The new- 
born pair of simple eigenvalues Xf corresponds to the forward and backward traveling waves, which propagate 
along the circumferential direction. Viewed from the stationary frame, the frequency of the forward traveling 
wave appears to increase and that of the backward traveling wave appears to decrease, as the spin increases. 
Double eigenvalues thus originate again at non-zero angular velocities, forming the nodes of the spectral mesh 
[lOj of the crossed frequency curves in the plane 'frequency' versus 'angular velocity'. The spectral meshes 
generically appear in equivariant dynamical systems fll and are characteristic for such rotating symmetric 
continua as circular strings, discs, rings and cylindrical and hemispherical shells, vortex rings, and a spherically 
symmetric a^-dynamo of magnetohydrodynamics [6l llOj. 

At the angular velocity flf — uj^/s the frequency of the sth backward traveling wave vanishes to zero 
(A J = A^ = 0), so that the wave remains stationary in the non- rotating frame. The lowest one of such 
velocities, ftcr, is called critical 12J. When the speed of rotation exceeds the critical speed, some backward 
waves, corresponding to the eigenvalues AJ , travel slower than the disc rotation speed and appear to be traveling 
forward (reflected waves). The effective energy of the reflected wave is negative and that of the forward and 
backward traveling waves is positive [13|. Therefore, in the subcritical speed region \Q\ < Qcr all the crossings 
of the frequency curves correspond to the forward and backward modes of the same signature, while in the 
supercritical speed region > flcr there exist crossings that are formed by the reflected and forward/backward 
modes of opposite signature. As a consequence of Krein's theorem |131 114| . under Hamiltonian perturbations 
the crossings in the subcritical region veer away into avoided crossings (stability), while in the supercritical 
region the crossings with the mixed signature turn into the rings of complex eigenvalues — bubbles of instability 
[13] — leading to flutter known also as the 'mass and stiffness instabilities' [6]. 

A supercritical flutter frequently occurs in the high speed applications like circular saws and computer 
storage devices, while the subcritical flutter — either desirable as a source of instability at low speeds as in the 
case of musical instruments like the singing wine glass and a glass harmonica or undesirable as in the case of the 
squealing disc- and drum brakes [3 [6l I12[ I15j — is an elusive phenomenon. The latter property is characteristic 
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for the regions with the definite Krein signature and appears, for example, as the subtleties in the tuning the 
helical turbulence parameter for excitation of the oscillatory MHD dynamo 

Being prohibited by Krein's theorem for the Hamiltonian systems, subcritical flutter can occur, however, 
due to non-Hamiltonian, i.e. dissipative and non-conservative, perturbations [6l[T7]. Since even the simplest 
codimension-1 Hamiltonian Hopf bifurcation can be viewed as a singular limit of the codimension-3 dissipative 
resonant 1 : 1 normal form and the essential singularity in which these two cases meet is topologically equivalent 
to Whitney's umbrella jTH [19l |20] , substantially complicating the stability analysis [211 IHl |23] , one can expect 
that exactly the singularities, associated with the double eigenvalues of the spectral mesh, make the induction 
of subcritical flutter in system ([1]) by non-Hamiltonian perturbations a challenging problem. 

In this Letter using perturbation theory of multiple eigenvalues we show that even if the eigenvalue branches 
in the subcritical region are well separated at the avoided crossings, created by the stiffness variation, they can 
be forced to bend with the origination of arcs of complex eigenvalues with positive real parts by the indefinite 
damping, which comes to the equations of motion, e.g., from the negative friction- velocity gradient [15112411^ . 
In the space of the gyroscopic, damping, and stiffness parameters the zones of the subcritical dissipation-induced 
flutter turn out to be cones with the apexes at the points, corresponding to the nodes of the spectral mesh. 
We show that the orientation of the cone is substantially determined by the structure of the damping matrix, 
which can be chosen in such a way that the system is unstable for significant magnitudes of the stiffness 
matrix detuning. The conical zones of the subcritical flutter bifurcate into couples of Whitney's umbrellas 
when non-conservative positional forces (originated from the follower forces or from the moment generated by 
variation in friction forces [12]) are additionally taken into account. We describe in detail this process related 
to the unfolding of the conical eigenvalue surface near a diabolical point to the "double coffee filter" singularity 
with two exceptional points under a complex perturbation of a real symmetric matrix in the problems of wave 
propagation in chiral and dissipative media [71 [H |26l [27l |28] . 

2 Conical zones of the dissipation-induced subcritical flutter 

The models of disc brakes, where subcritical flutter was detected by numerical approaches, include both the 
case when the point-wise or distributed friction pads are rotated around a stationary disc, affecting a point or 
a sector of it, and when the disc rotates past the stationary friction pads, see [SI I12j and references therein. 
Linearization and discretization of the latter class of the models frequently results in the equation ([T]) perturbed 
by the matrices D = (dy) — T)^ , K = (hj) = K-'", and N ~ (riij) = — N-^ of dissipative, potential, and 
non-conservative positional forces, respectively: 

X + (20G + (5D)i + (P + n^G^ +kK + i^N)x = 0, (5) 

where the parameters 5, k, and v control the magnitudes of the perturbations. The matrices D, K, and N can 
be assumed to be functions of Q. 
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Figure 1: Numerically calculated eigenvalue branches (thick lines) and their approximations ([5]) (dashed lines) 
for uji = 1, fell = 1, k22 = 2, fci2 = 1, dii = -1, c?22 = 2, di2 = 0, J/ = 0, K = 0.2, and (a) 5^0, (b)-(d) 5 = 0.3. 

In the rotating frame the load appears to be moving periodically in the circumferential direction. The 
transformation x = Az := exp(— riGi)z, see [IHl 133 , yields an equivalent to © potential system with the 
time-periodic perturbation of the coefficients 

z + (5D(<)z + (P - 5nf){t)G{t) + K'k{t) + iyN(t))z = 0, (6) 

where T>{t) = A^^DA, G{t) ^ A^^GA = G, K(t) = A^^KA, N(t) = A^^NA. For 2n = 2 degrees of 
freedom the matrix N(t) = N, while the periodic stiffness and damping matrices have a simple explicit form 

2K(i) = diag(trK,trK) + (K + JKJ)cos(2m) + (JK-KJ)sin(2m), 

2D(i) = diag(trD,trD) + (D + JDJ)cos(2m)-|-(JD-DJ)sin(2m). (7) 

Note that in recent Letters [30l [3T] , Shapiro studied only a particular case of system ([6]) , ([7]) with 2n = 2 degrees 
of freedom, D = diag (1, 1), K = diag (—1, 1), = 3 , lo\ — a, n — 2q, 5 — 2h, and i' — 2hn, which is, however, 
used much in modeling strong focusing in electron and ion optics through a twisted series of static quadrupole 
magnetic or electric lenses or helical quadrupole channels as well as in modeling the monochromatic light of two 
orthogonal polarizations propagating along the axis of helically twisted anisotropic media like cholesteric liquid 
crystals, see [301 ISI] and references therein. 

In the absence of dissipative {S — 0) and non-conservative positional {v — 0) terms, (jH) describes coupled 
parametric oscillators of Mathieu-type with the periodic in time potential possessing parametric resonance in the 
supercritcal range > flcr- Inclusion of parametrically excited damping and non-conservative terms makes 
the equation ^ a less traditional periodic system demonstrating parametric resonance in both the subcritical 
and supercritical regions [12]. 

Nevertheless, the equivalence of the two dual descriptions enables us to reduce the investigation of the 
periodic system ([B]) to a considerably simpler study of the stability of the autonomous system (O, cf. I30j . 
In particular, subcritical parametric resonance domains of equation ([6]) correspond to the regions of subcritical 
flutter of the system ([5]). 
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Figure 2: Domains of the subcritical flutter instability for indefinite damping matrix in the absence of the 
non-conservative positional forces when (a) A > 0, ih) A — 0, and (c) A < 0. 



Owing to its relative simplicity the case of two degrees of freedom (n = 1) allows for the detailed stability 
analysis. For n — \ the spectrum of the unperturbed operator Lo(f2) consists of four branches (U, forming 
the spectral mesh in the (£7, ImA)-plane. In the subcritical region < flcr = the branches cross at the 
points (0, ittJi). Assuming without loss in generality N = J, we consider a general perturbation of the gyroscopic 
system Lo(ri) + AL(fi). The size of the perturbation AL(ri) = JAD + kK + z/N ~ e is small, where e = |jAL(0)|| 
is the Frobcnius norm of the perturbation at Q = 0. For small Q and e perturbation of the double semi-simple 
eigenvalue A = icui with two orthogonal eigenvectors Ui — (0, (2a;i)^^/^) and U2 = ((2a;i)^^/^, 0) is described 
by the formulas [6] 



c -I- Rec 



IniA = LUi 



Pi + P2 
4wi 



K± 



Rec 



Rec 



Imc = 



6k- 



2trKD trKtrD 



8cji 



(8) 

(9) 



4 / V 4cji 

where /ii.2 and pi^2 are the eigenvalues of the matrices D and K, respectively. 

Since the eigenvalues at the crossings in the subcritical range have the same Krein signature, they veer away 
under potential perturbation kK, destroying the rotational symmetry of the body, independently on whether 
the perturbing matrix K is definite or indefinite [6] . The detuning of the stiffness matrix P + kK does not shift 
the eigenvalues from the imaginary axis, preserving the marginal stability. Near the node of the spectral mesh 
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(ri = 0, ImA — uJi) the veering is described by the hyperbola [3] 



Pi + P2 A o2 Pi- P2 \ 2 
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ImA - c^i - ^li-^'* = rV^ : ReA-0, (10) 

shown by the dashed hnes in Fig. [Ija). 

Now we demonstarte that even if the eigenvalues of the rotationally symmetric gyroscopic system ([T]) are 
separated in the subcritical region by the symmetry- breaking variation of the stiffness matrix kK, the inclusion 
of dissipation 6T) can cause flutter instability when D is indefinite, see Fig. [Hb)-(d). Indeed, with = in 
equation ([8]), the condition RcA < yields linear approximation to the domain of asymptotic stability in the 
(fi, K, 5) - space 

(5trD>0, K^A + Vl?{2u)iiYT)f > -det'D{ujitvT>f5^, (11) 



where 

A:=detD(pi-p2)'+(A:i2(d22-dii)-di2(fc22-A:ii))^- ^^^^^^^ ^"7^^'^^ , (12) 



,2 ((trD)2-16/32)(pi_p2)2 



2trKD-trKtrD , , 

/3o 77 ^ • (13) 

4(pi - P2) 

For the definite damping matrices we have A > because det D > 0. Consequently, for 5T) > the conditions 
(fTTj) are always fulfilled in agreement with the Thomson- Tait-Chetaev theorem [221 US]- Since detD < in case 
of indefinite damping, the inequalities ([TT|) indicate that the flutter instability domain for A > is inside the 
conical surface extended along the (5-axis and the stability domain is adjacent to the cone's skirt selected by the 
condition (5trD > 0, see Fig. [2][a). The conical domain is stretched along the K-axis when A tends to zero and 
it is transformed into a dihedral angle when A = 0, as shown in Fig.[2{b). With the further decrease in A the 
dihedral angle is again wrapped into the conical surface which is then extended along the fi-axis, Fig.[2)Jc). The 
domain of asymptotic stability is inside the half of the cone selected by the inequality ^trD > 0. The eigenvalue 
branches with their approximations ([S]) , illustrating the mechanism of the dissipation- induced subcritical flutter 
in the presence of imperfections kK, are shown in Fig. [Ijb)-(d) for A < 0. Note that the threshold ^4 = 0, 
separating the indefinite damping matrices, coincides with that found first in [22| from the criteria of Routh 
and Hurwitz for a general two-dimensional non-conservative gyroscopic system with dissipation. 

In case when n > 1, the mutual orientation of the conical stability boundaries is expected to be determined 
also by the Krein signature of the eigenvalues involved into the corresponding crossings at the nodes of the 
spectral mesh, which is substantially different in the subcritical and in the supercritical regions, similarly to the 
Krein space induced geometry of the resonance tongues of the spherically-symmetric MHD a^-dynamo 16J. 

In the plane (fi, k) for a fixed 6 > the instability domain has, respectively, the form of an ellipse, Fig.l^I^a), 
a stripe. Fig. [HJb), or a region contained between the branches of a hyperbola, Fig. HJc). The latter picture 
shows that a widely known in the engineering practice approach to the squeal suppression by the reduction of 
the rotational symmetry of the rotor is not efficient in the presence of indefinite damping, which originates from 
the brake pads with the negative friction- velocity gradient [HI [24l ES] . And daily experience testifies, that an 
imperfect wine glass and a perforated brake disc still experience friction-induced instabilities and emit sound. 
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Figure 3: Unfolding the conical boundary of subcritical flutter (dashed lines) by the perturbation i^N into a 
couple of Whitney's umbrellas with their level curves depicted for 5 > when (a),(b) ^ > and (c),(d) ^ < . 

Bifurcation of the stability diagrams with the change of the entries of the matrices D and K can also explain 
unsatisfactory reproducibility of experiments with disc brakes [T2]- Indeed, some parameters like rotational 
speed and pressure on the brake pads can be regulated precisely, while the topography of the pads' surface as 
well as the material properties of the pads undergo uncontrollable changes from one run of the experiment to 
the other owing to the heating, cooling, and wear |12j . As a consequence, the very same values of stiffness n 
and angular velocity 51, which yield flutter instability and squeal for < 0, Fig.[5]Jc), make the brake quiet for 
A> Q owing to the qualitative change of the stability diagram with the change in the structure of the matrices 
D and K, Fig.[2^a). This effect becomes even more pronounced in the presence of non-conservative positional 
forces, because, as we demonstrate in the next section, the conical stability boundaries are not structurally 
stable under small perturbations i^N. 



3 Unfolding the conical zones into couples of Whitney's umbrellas 
in the presence of non-conservative positional forces 

In the presence of non-conservative positional forces (j/ 7^ 0), stability conditions (jlip . following from equations 
([HI) under the requirement ReA < 0, are modified into (5trD > and i? > 0, where 

B := (2n((52w^trD2 - 4j/2) + 5(2trKD - trKtrD)Ki/)2 + ShrB^iAS^uf - iy^{pi ~ p2f)K^ (14) 

- Shr'D^iS^LoltT'D^ ~ - S'^ujf det D). 

In the assumption that u — the inequality i? > is reduced to the second one of (fTT|) : for k = it is 
transformed into fl^ < il^,., where [6] 

_ (5trD / i^^ ~ ujf 5^ detB 
^^-■-^y-^,2„^2^2(trD/2)2- ^^^> 

For the fixed parameters S and i' the equation B = generically describes either an ellipse or a hyperbola 
in the (17, K)-plane, Fig. [3jb),(d). For A> and sufficiently big S, the domain of subcritical flutter is inside 
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the ellipse B ~ 0, and the domain of asymptotic stability is outside the ellipse similarly to the system without 
non-conservative positional forces, whose stability diagram is shown in Fig.l^^a). With the decrease in S for the 
fixed v the ellipse is rotated around the origin in the (il, K)-plane and simultaneously it is stretched along one 
of its main axes. At the threshold AS'^cof = v^{pi — p-if' the ellipse is transformed into two parallel lines, which 
with the further decrease in 8 are bent into two branches of hyperbola i? = 0. Then, the flutter instability 
domain lies between the two branches, Fig.[31[b). When 5 tends to zero, the hyperbolic branches are shnmk to 
the miion of the intervals (— cxd, — kq] UI'^o, oo) of the K-axis, where the critical value 

"^v , ^ 

Ko := 16 

Pi - Pi 

follows from the conditions Rec — and Imc = implying existence of the double eigenvalues at the points 
(17 = 0, K = ±Ko, (5 = 0), which are shown as open circles in Fig. [31 

Remarkably, the value (jl6p of kq obtained from the perturbation formulas ([8]) coincides with the exact one 
following from the characteristic equation of the system ([5]) in the assumption that (5 = and il = 

A* + (2u:l + KtrK)A2 + det K + Kw^trK + = (17) 

Substituting k = kq into (jl7p yields the frequency of the corresponding double eigenvalues ±iwo 



UJo-=\h^t + l' =UJi + - ho(j^), (18) 

Pi - P2 2UJI pi - P2 



where the first-order approximation with respect to v follows also from the formulas ([8|) and p6p . 
The double eigenvalue iujo has the eigenvector uq and the associated vector ui of the Jordan chain 

\^ 2fci2 +P1-P2 J " \0 J 

which are the solutions to the equations 

(-w^I + P + KoK)uo = 0, (-w^I + P + KoK)ui = 2iwoUo. (20) 

Therefore, in the parameter space the coordinates (0, ±ko,0) correspond to exceptional points (EPs), at which 
there exist double eigenvalues with the Jordan chain. Fig. [3l 

In the vicinity of the EPs the expression B — for the stability boundary yields 



n = ^^^^^l^^^^^^^^^s + oiS). (21) 

Thus, in the {Q, (5)-planc the domain of asymptotic stability is bounded in the first-order approximation by the 
two straight lines ((2T|) . When k goes to ±kq, the slopes of the both lines [3 :— Vl/6 tend to the value (3 — ±/3o, 
where /3o is defined in p3p . Extracting k from the equation (j2ip we find another representation for the stability 
boundary near exceptional points 



4/3/3o±trD^/32-/32 + (i^)^ 



1 + : 



trD 



+ o{il3Tl3of), (22) 



8 



which has a canonical for the Whitney's umbrella singularity form Z = X'^/Y'^ ^TE[ [TOl [22] , 

Therefore, we explicitly demonstrated that the conical boundary of the domain of subcritical flutter for 
A > is structurally unstable under the perturbation i/N. With the increase of v the cone opens up and 
simultaneously the plane 6 = foliates into two sheets intersecting along the branch cuts (±cxd, ±ko] on the 
K-axis, which are shown as thick lines in Fig.[3l The new surface has a couple of Whitney's umbrella singularities 
at the exceptional points (0, ±ko,0). The domain of asymptotic stability, which for v = was adjacent to the 
conical domain of subcritical flutter is now wrapped into the pockets of the two Whitney's umbrellas, selected 
by the inequality (5trD > 0. With the increase in 6 the stability boundary gradually tends to the conical surface 
with the flutter instability inside it. Fig. ^a). 

Inclusion of the non-conservative forces qualitatively changes the stability diagram in the (17, K)-plane trans- 
forming the elliptic flutter domain of Fig.[2lja) to a larger one located between the hyperbolic branches, Fig.[3l^b). 
For 1^ = and (5 = the {fl, K)-plane is stable, while for Q and (5^0 the stability domain dramatically 
shrinks to the branch cuts (±oo, ±ko]. Consequently, under small perturbation vN a point in the (fi, K)-plane, 
which was in the stability domain for v = can suddenly find itself in the instability region when v ^ 0, 
similarly to the scenario described in the previous section. 

For A < the conical stability boundary of Fig. ^c) unfolds into two surfaces with the Whitney's umbrella 
singularities at the exceptional points (0, ±ko, 0) as shown in Fig.[3]Jc). The local approximations to the surfaces 
near the singularities are given by the same equation (j22p . where /3o has a value different from the case when 
A > 0. For i> = and 6^0 the stability domain in the {Q, K)-plane is inside of the two hyperbolic regions 
extended along the fi-axis, as shown in Fig.jJ^c). When ^ with the decrease of 6 the stability domain rotates 
around the origin until it is completely reoriented and shrunk into the branch cuts (±oo, ±ko] extended along 
the K-axis, Fig. [31(d). Due to such a reorientation one can again observe sudden stabilization/destabilization at 
the very same values of fl and n when v is slightly changed. Similar discontinuous change of planar stability 
diagrams was observed in [THl [301 [31] for parametric oscillators coupled through the gyroscopic and damping 
forces. 

Therefore, the fundamental process of birth of a couple of exceptional points from a diabolical point (whose 
location is marked by the asterisk in Fig. [3]) , which has important consequences for the wave propagation in 
chiral and dissipative media [71 [51 [^71 [^ , substantially determines the stability of waves propagating in rotating 
continua with frictional contact. 

Conclusions 

The singing wine glass and the squealing brake are phenomena of the acoustics of friction, which everyone 
encounters with almost every day. Despite of a seeming simplicity, their mathematical modeling is not easy, 
and experiments with them are not satisfactorily reproducible. Several fundamental reasons are responsible for 
this. First, the rotational symmetry of the continua makes the spectral mesh with the double eigenvalues at the 
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nodes generic for the unloaded bodies of revolution. Then, definite Krein signature at the eigenvalue crossings in 
the subcritical rotation speed range requires dissipative and non-conservative perturbations to produce complex 
eigenvalues with the positive real parts. The geometry of the domain of the dissipation-induced subcritical flutter 
is complicated by the singularities related to multiple eigenvalues. The changes in the damping and stiffness 
matrices, which are unavoidable from one run of the experiment to another one due to wear of the friction 
pads, provoke bifurcation of the flutter domain, which explains why experimental results are not reproduced 
even after the accurate restoration of other well-tunable parameters. Using perturbation theory of multiple 
eigenvalues we described typical bifurcations of the subcritical flutter domain related to the birth of a couple 
of exceptional points from a diabolical point due to non-Hamiltonian perturbation of a Hamiltonian system. 
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